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Abstract 

This is a proof of an asymptotic formula which describes exponen- 
tially small splitting of separatrices in a generic analytic family of area- 
preserving maps near a Hamiltonian saddle-centre bifurcation. As a 
particular case and in combination with an earlier work on a Stokes 
constant for the Henon map (Gelfreich, Sauzin (2001)), it implies ex- 
ponentially small transversality of separatrices in the area-preserving 
Henon family when the multiplicator of the fixed point is close to one. 
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There has been a substantial progress in the study of exponentially small 
splitting of separatrices (see for example [H [21 EEH1 EB 122])- In most cases 
exponentially small asymptotic expansions are very sensitive to the form of 
the system. In the case of local bifurcations with an integrable normal form 
some generically valid formulae have been obtained p^8l fTTl [T2] . In degenerate 
cases the splitting of separatrices may be studied by the Melnikov method 
[21 H] but application of this method to exponentially small phenomena add 
substantial restrictions onto systems under consideration (see e.g. survey 
[T3]). Since the pioning work by Lazutkin, it is accepted that in the case when 
Melnikov method fails the asymptotic formulae contain a special pre-factor 
which comes from a parameter-independent problem and can be interpreted 
as a Stokes constant (see e.g. [IS])- For the class of system considered in the 
present paper, that problem is considered separately in [T3] . 

We note that transversality of separatrices plays an important role in 
some applications ([El [7]). 



1 Introduction 

Let f e be an analytic family of area-preserving maps defined in a neighbour- 
hood of the origin. We assume that when e = the origin is a parabolic fixed 
point of the map. This means that /o(0) = and one is a double eigenvalue 
of the Jacobian Df (Q). These eigenvalues are called multipliers of the fixed 
point. In the generic case the Jacobian is not diagonalisable. Then we may 
assume that 

£/o(0) = ( J \ ) , (1.1) 

as this form can be achieved by a linear area-preserving change of variables. 
Since the multipliers are equal to one the fixed point is not structurally stable 
and can be destroyed by a small perturbation of /q. We will consider the 
non-degenerate bifurcation. This means that two leading coefficients in a 
normal form of f £ do not vanish: one of them is a leading non-linear term of 
fo and the second one guarantees that f E is a "generic unfolding" of fa. 

Let us state these two assumptions in a more precise way. We will see 
that the bifurcation can be described by an integrable normal form. We will 
provide a detailed description of the normal form later but at the moment 
we note that the map f £ can be approximated by a time one map of an 
autonomous Hamiltonian flow. The leading order of the Hamiltonian has 
the form 

2 3 

h& = y — + — -bxe. (1.2) 
6 2 3 y ' 
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Figure 1: Bifurcation of the normal form: a small separatrix loop is created 
after a Hamiltonian saddle-centre bifurcation 

In this paper we assume 

a,b>0. (1.3) 

The positivity assumption is not more restrictive than ah ^ 0. Indeed, the 
sign of a and b can be changed by substitutions x \— > — x and e \— > — e. Using 
assumption (11.31) we ensure that the phenomena we study in this paper are 
on the side of positive e. 

Near the bifurcation the map f e is rather accurately approximated by the 
time one map <&\ generated by the Hamiltonian (jl.2p . When e changes sign 
an equilibrium of h% bifurcates following the scenario sketched in Figure HJ 
For e > 0, separatrices of the saddle point form a loop around the elliptic 
fixed point. 

Under these assumptions f e with e > has two fixed points which collide 
at the origin when e = and then disappear. One of the fixed points is 
elliptic (centre) and another one is hyperbolic (saddle). For this reason the 
bifurcation is sometimes called a Hamiltonian saddle-centre bifurcation. 

The bifurcation of the map is very similar to the bifurcation of the lead- 
ing order normal form (II. 2J) . Unlike the normal form, we do not expect a 
generic map to be integrable and therefore the separatrices of the map may 
split. We will derive an asymptotic formula which describes the splitting of 
the separatrices. This formula implies that in a generic family of maps the 
separatrices intersect transversally for all sufficiently small e. 

We note that the area preserving Henon map satisfies assumptions of our 
theorem. 

The leading order of the asymptotic expansion was described in [TT| [T2"] . 
The goal of this paper is to present a complete proof of a refined version of 
this asymptotic formula. The proof is based on the programme of [TO] and 
is a development of the original idea proposed by Lazutkin in |17j . 

Let us state the main result of the paper. First we note that many of the 
series involved in our proof involve powers of e 1 ^. It is therefore convenient 
to introduce a new small parameter 5 = e 1//4 . We will occasionally use both 
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of them indicating the type of a power series expansion for the corresponding 
function. 

In this paper we will study separatrices of a family of hyperbolic fixed 
points: 

P<5 = fe(Ps) 

such that p = 0. Let \$ > 1 be the largest of the two multipliers of f £ at 
p<5. We will prove that p<s and Xs depend analytically on 5. Moreover, logA,s 
is of the same order of smallness as 5 itself. We will see that 

logA 5 = \^Aa~b5 + 0{5 3 ). (1.4) 

It is convenient to represents stable and unstable separatrices W /± (p,5) 
using solutions of the following finite-difference equation (for a more detailed 
discussion see the review [13]) 

V> ± (t + logA (5 ) = / e o^(t). (1.5) 

We have suppressed explicit dependence of ^ on 5 to shorten the notation. 
An application of f £ increases the parameter t by logA^ > 0. Therefore the 
conditions 

lim ip~(t) = ps and lim ip + (t) = p$ 

imply that ip~ represents the unstable separatrix and ip + represents the stable 
one. It is easy to see that these conditions do not defined the functions V ;± 
uniquely. 

The most convenient parameterization is obtained using analytic lin- 
earizations of restrictions of f E onto W~(ps) and W + (p$) respectively. These 
parametrisations can be represented in the form ip~(t) = ^f~(e t ) and ip + {t) = 
\l/ + (e _ *) where ^/ ± (z) are analytic at and \E f± (0) = p,5. This choice reduces 
the freedom in ip^ to a translation 1 1— > t+to- Our proof of their existence does 
not use the linearization but relies on a contraction mapping argument which 
has the advantage of providing an accurate approximation for the functions. 

From the technical point of view a substantial part of the present paper 
is dedicated to a detailed study of the analytic continuation for ip ± . 

Let 7 e W~(p$) fl W + (ps) be an homoclinic point. Then 

7 = ^+(t s ) = 

for some t u and t s . The vectors ip~(t u ) and ip + (t s ) are tangent respectively 
to ^"(p^) and W + (ps) at 7. The area of the parallelogram generated by 
these tangent vectors is called the homoclinic invariant 0/7: 

w( 7 ) = ft (<Ht s ),^ + (t u )) , 
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where f2 stands for the standard area form. The homoclinic invariant has 
some advantages over alternative measures for the separatrices splitting (see 
the review [13] for a discussion). 

Theorem 1.1 (Main Theorem). There are constants £ M. such that the 
homoclinic invariant of each of the two primary homoclinic orbits has the 
form 

u^±-^e-T&Ta k 5 2k . (1.6) 

This theorem implies that for a generic family the separatrices intersect 
transversally. Indeed, we note that in (II. 6p the pre-factor is exponentially 
small compared to 5 due to the smallness of log A<5 estimated in (11.4j) . Never- 
theless for small 5 > the theorem implies transversality of the homoclinic 
points provided a ^ 0. We note that a describes the splitting of complex 
parabolic invariant manifolds for function / and does not vanish generically 
(on an open dense subset). In [15J it was proved that a 7^ for the Henon 
map. Finally, numerical evaluation of ao for a given fo is a relatively easy 
task. 

The difficulty of this theorem is due to the exponential smallness of the 
separatrices splitting which is to be detected on the background of much 
larger effects. 

The proof consists of the following main steps: 

1. We show that f £ can be formally represented as a time-one map of a 
Hamiltonian flow and provide a normal form theory for the flow. 

2. We develop a general theory of close-to-identity maps which provides 
sufficient condition for the existence of saddle points and their separa- 
trices. Moreover, we prove that the separatrices of the map can be quite 
accurately approximated by separatrices of a flow which approximately 
interpolates the map. 

3. We show that In the complexified time domain the separatrix of the 
interpolating flow is close to the separatrix of f £ in a set which has a 
non-empty intersection with a ^-neighbourhood of ±z7r. 

4. We show that in the latter region an approximation based on the sepa- 
ratrices of the parabolic fixed point of fo provides an even more accurate 
approximation. This approximation distinguishes between the stable 
and unstable separatrices of f £ . 
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5. We show that flow-box time coordinates can be introduced. 

6. Flow-box coordinates are used to get an asymptotic expansion for the 
splitting of separatrices. 

The rest of the paper is dedicated to proving the main theorem and is 
structured in the following way. Section [2] contains a detailed overview of 
the proof leaving technical estimates to corresponding sections. Section [3] 
describes the construction of the normal form for the bifurcation up an ar- 
bitrary order. Section develops a general theory of close to identity maps, 
which can be of independent interest since it covers a much wider class of 
maps than that studied in other parts of this paper. In Section H] we provide 
a detailed description of a formal solution of the separatrix equation (11.51) . 

We will prove that f e can be analytically interpolated by an autonomous 
analytic flow and establish properties for the corresponding energy-time co- 
ordinates. Note that this property does not imply integrability of f £ since the 
domain in question is not invariant. Several sections are included to provide 
details on approximation of ip~{t) on various subsets of the complex plane. 

2 Overview of the proof 

In this section we will describe the main steps of the proof postponing proofs 
of technical statements to subsequent sections. Before proceeding to the 
proof we need to introduce some definitions. 

In this paper we will use power series in x, y and e. The arguments are 
easier when terms of a similar magnitude are grouped together. We consider 
x to be of order 2, y to be of order 3 and e of order 4. We say that a 
polynomial Q p is quasi-homogeneous of order p if it satisfies the identity 

Q p (X 2 x,X 3 y,X 4 e) = X p Q p (x,y,e). (2.1) 

Any power series can be represented as a sum of quasi-homogeneous polyno- 
mials. 

In order to shorten the notation, we define 

d n :=0[(max{\x\^M 1/a ,^}) n ) 

and write O n , n +i for a vector with components (O n ,O n+ i). It is easy to see 
that Q p = O p . 
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Figure 2: Overview of the proof: Domains of approximation 
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2.1 Standard scaling and limit flow 

The change of the variables 

x = 5 2 X, y = 5 3 Y, e = 6 4 (2.2) 

is called the standard scaling. In the new variables the map f e takes the form 
of a close to identity map 

F s = id + 5G 5 . (2.3) 

We note that Go is called the limit flow. Its time-5 map & 6 Go = id+5Go+0(6 2 ) 
and therefore is 0(5 2 ) close to F$. 

Taking into account that /o(0) = and the equality fll.ip we obtain 
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Y 
-aX 2 



for some constants a, b. The same constants enter the non- degeneracy as- 
sumption (11.31) . Indeed, the limit flow is Hamiltonian with the Hamilton 
function 

2 3 

which coincides, after reverting the standard scaling and scaling time, with 
(11.21) . The dynamics of the limit flow is defined by the system of Hamiltonian 
equations: 

X = Y, Y = -aX 2 + b. 



The limit flow has a saddle equilibrium at (— y ^ ; 0). At this point the lin- 
earised vector-field has a positive eigenvalue fi = \/ Aab. The separatrix of 
the limit flow can be found explicitly by integrating the Hamiltonian equa- 
tions: 



, ^ lb 3Vab . . r- sinh^ 



a ' cosh 2 if 1 uv ' ^ u C osh 3 ^ ' 

A substitution shows that these functions satisfy the Hamiltonian system. 
Obviously they converge to the saddle equilibrium as t — ► ±oo. 

It is convenient to scale time and define the separatrix solution of the 
limit flow by 

<p (t) = (X ,Y )(t/fx ). (2.4) 
This function satisfies the differential equation 

0o = 1 JVi? (^o) • 
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We will consider the separatrix solution for complex values of t. We note that 
ip is analytic for all complex t except for second order poles at iir(2k + 1) 
with k G Z. In particular this solution is analytic in the strip |Im(t)| < 7r 
and has poles at the points = ±7ii on its boundary. These singularities 
play the central role in our analysis. 

2.2 Formal interpolation and Formal separatrix 

We use h™ = £™=g h p (x, y, e) to denote a partial sum of a formal series 



where h p (x,y,e) are quasi- homogeneous polynomials of order p. In Sec- 
tion 13.21 we will prove that there exists a unique formal series such that for 
each n e~N 



where Q\ n is the time one map with hamiltonian /i™. The leading order of the 
series is given by ( 11.21) . It is a quasi- homogeneous polynomial of order 6 and 
coincides, up to scaling of the space and time variables, with the Hamiltonian 
of the limit flow. 

We say that the series h e formally interpolates the map f e . We should 
not expect the series h £ to converge. 

Theorem 2.1 (Formal separatrix). Ifh £ is a formal Hamiltonian of the form 
(12.51) and its leading order satisfies the non- degeneracy condition (11.31) . then 
there is a unique formal series 




(2.5) 



fe — ®h™ + ^n+3,n+4 , 



(2.6) 



■2k+l 



such that the equation 



/i(5)X = JV/i e (X) 



(2.7) 



has a formal solution of the form 



p>i p>i 



X(t,e) 



(2.8) 



p>l p>2 



where X*, Yp,Xp, Y p 2 are polynomials of order p in i] = cosh 



-2 t 
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The proof of this theorem is placed in Section HI We remind that 5 = e l / A . 

We say that X is a formal separatrix. The series h £ may diverge and con- 
sequently there is no reason to expect that the formal separatrix converges. 
On the other hand we will see that its partial sums provide a rather accurate 
approximation for the separatrices of the map f £ . 

Let us consider the domain 

T = { t G C : Re(t) < n and <p (t - s) G V, Vs > } (2.9) 

where r\ > is an arbitrary constant and T> is a bounded domain. Later we 
will assume it to be a sufficiently large ball centred around the origin. Of 
course, some constants in future estimates will depend on the choice of r\ 
and V. 

Lemma 2.2. For any n G N, the Hamiltonian h$ has a saddle equilibrium 
with a Lyapunov exponent fi Uy s > in a neighbourhood of the origin. More- 
over, there exists a separatrix solution tpg of the equation 



such that uniformly for t G % 

where X n denotes the sum of the first n orders in 6 of the formal series 



Proof. The proof is straightforward: this lemma is a statement about poly- 
nomial vector fields and after the standard scaling the saddle of the equation 
is non-degenerate. □ 

Let X be the formal separatrix given by (12. 8p . We will study it close to 
the singularity by substituting t = m + r log \$ and expanding Xl' 2 , Y k ' into 
their Laurent series. We also re-expand logA<5 = //(£). The result is a formal 
series of the form 

V I l2mJm-2 \ ^ ^ -k \^ r2mJm-3 \ ^ „ —k \ /<-> i i \ 

\m>0 k>0 m>0 k>0 J 

where x m k and y m k are real coefficients. See Section 14.31 for a derivation of 
(12. lip . It is convenient to introduce formal power series in r by setting 

= [r 2m - 2 J2 x m kr~ k ,T 2m - 3 J2ymkr~ k J (2.12) 

\ k>0 k>0 J 
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for m > 0. Then 

X=^5 2m ^ m (r). (2.13) 

m>0 

We will use these formal series in the complex matching procedure described 
in Section [2761 



2.3 Separatrices for close-to-identity maps 

Let us consider a family of close to identity maps (12.31) . For the purpose of 
this section it is not necessary to assume that the map is obtained as a result 
of the standard scaling. 

The implicit function theorem implies that if the limit flow Go has a non- 
degenerate equilibrium then F$ has an analytic family of fixed points which 
tend to the equilibrium when 5^0. If the equilibrium is a saddle then the 
fixed point of F$ is also a saddle (see Section [5711 for formal statements and 
proofs) . 

The separatrix of F$ is close to the separatrix of the limit flow (see [9]). 
We will use a more accurate approximation provided by a separatrix of an 
(approximately) interpolating flow. Let us state it more formally. 

Let T> C C 2 be a bounded domain and assume the origin is inside T>. 

Theorem 2.3 (Approximation theorem). Let F$ be an analytic family of 
area preserving maps such that the limit flow has a saddle equilibrium at the 
origin. Let H$ be a Hamiltonian function such that 

F s = + 0{5 n+1 ) (2.14) 

6 

on T> and be a separatrix solution of the equation 
such that 

<%(t) = <po(t)+0(S) 

for all t G T . Then there exists a parametrisation \l/~ of the local unstable 
separatrix of the map Fs which satisfies the equation 

*-(t + log\s) = F s (*-(t)) (2.15) 

and 

*-(t)=$(t)+0(5 n ) (2.16) 

uniformly on the set T . 
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The proof of this theorem consists of two steps. First the theorem is 
proved for the local separatrix which corresponds to a half-plane Re t < —r 
as described in Section 15.31 Then an extension lemma (Lemma 15.31) stated 
and proved in Section 15.41 is applied. 

For completeness, we also provide a proof of the existence for the Hamil- 
tonian in Section 15.21 

2.4 Parametrisation the separatrices 

We note that the existence of a parameterisation which satisfies equation 
(11.51) follows from the general theory [13] . The parametrization is not unique 
and is defined up to a translation 1 1— > t + t (e). The next theorem states that 
there is a parameterisation which is close to the separatrix of the interpolating 
flow. 

Let be defined by Lemma 12.21 

Theorem 2.4 (Existence and Local approximation). Equation ( II. 5\) has a 

solution 

tf;-(t) = ^(t)+O n+2 , n+3 (S) (2.17) 

uniformly for t G 7q. 

Proof. The theory of the previous section applies to the map Fg obtained 
after the standard scaling of the map f £ . Let T> be a large ball of a fixed 
radius. Since the map f e is defined in an e independent neighbourhood of 
the origin, the domain of the scaled map F$ contains the ball V provided e 
is sufficiently small. 

The limit flow separatrix ip is given explicitly by (12. 4p and has poles at 
in (2k + 1) for each k e Z. Therefore the set Tq takes the shape similar to 
the one shown on Figure [3j it is a half plane Re t < ri without "shadows" of 
small disks centred around the singular points of (po. 

The interpolating Hamiltonian is obtained from (12.51) after the change of 
variables (12.21) . Although the change of variables is not symplectic it has a 
constant Jacobian, and therefore the Hamiltonian flow of is transformed 
into a Hamiltonian flow with the Hamiltonian function 

n 

H% = Y i 6 k h B+k (X,Y,l). (2.18) 
fc=i 

Indeed, under the standard scaling the symplectic form becomes dx A dy = 
5 5 dX AdY , to work with the standard symplectic form in the scaled variables 
we scale the Hamiltonian by 5~ 5 . The upper bound (12.141) . which is necessary 
to apply Theorems 12.31 follows from (12. 6ft and (12. 2\\ . 
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Im(t) 



Re(t) 



Figure 3: Domain T . 

To complete the proof we reverse the standard scaling which transforms 
the error term 0(5 n ) of equation (I2.16P into O n+2 n+3 ((5). □ 

2.5 Extension towards the singularity 

The theory presented in the previous subsection provides a rather accurate 
approximation for the invariant manifolds. However, it is not sufficient to 
distinguish between the stable and unstable separatrices of the map F$. In- 
stead of further improving the accuracy, we will study the separatrices for 
values of t closer to the singular points of the limit flow separatrix. In this 
region the splitting of the separatrices is larger and consequently easier to 
detect. This extension is sensitive to the form of the map and cannot be 
performed with the same level of generality as the estimates of Section 12.31 
So from now on we restrict our consideration to our original map f E . 

Let us fix two small positive constants c\ and /3, and consider a domain 
Ti{e) shown on Figure H] (left), which also depends on a constant c 2 > 0. The 
constant c 2 is to be chosen sufficiently large to ensure that (po, the separatrix 
of the limit flow written in the unsealed variables, does not leave a small 
e-independent ball centred around the origin. 

It is convenient to work with a time parameter centred at the singularity 
and to scale the time step to one, hence we let 

t — in . 
r=-— . 2.19 
log X s 

The domain Ti(e) expressed in the terms of r is shown on Figure H] (right). 
Theorem 2.5. Let ip~ be defined by Theorem \2.4\ There is c 2 > such that 

tfj-(T) = ^(r)+O n , n+1 (r- 1 ) (2.20) 
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Figure 4: An amplification of a neighbourhood of t = m shows the domain 
Ti{e) in which the unstable separatrix is still reasonably well approximated 
by the separatrix of the limit flow (left). The same domain on the plane of r 
(right). 



uniformly for all r G Ti(e) and e G (0, £o)- 

We note that t = m + r log Xs and we have overloaded notation by writing 
ip~ (r) instead of ip~ (in + r log Xs) = ip~ (t) . 

The estimate H2.2QI) is uniform which means that the constant in the term 
0„ iri+ i(r _1 ) is independent from both e and r. It is important because the 
size of 7i(e) increases as e — > and includes r up to the order of 5' 1 = e~ l l A . 
We note that in Theorem 12 .51 the error of the approximation is of order 0(S n ) 
on the left boundary of %_{e) but it is just of order 0(1) on the boundary 
near the central circle. 

Therefore the theorem suggests that the distance between the separatrix 
of the map and separatrix of the interpolating flow may gradually increase 
as the parameter t comes closer to in. 



2.6 Complex matching 

In this section we construct another approximation for ip~ which provides 
higher accuracy in the central part of T\{s). Let us start by looking for a 
formal solution 

^r) = Y J t 2k Mr) (2.21) 

k>0 

of the separatrix equation 

$(T+l) = fS(T)). (2.22) 

We assume that the new time r is related to the original t by (12.191) . We 
formally substitute the series (12.211) into the equation (I2.22p and expand the 
right hand side in powers of 5 2 . Collecting terms of equal order in S we get 
at the leading order 

Mr+l) = fo(Mr)), (2.23) 
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and the following system of equations for all other orders: 
Vi(r + 1) = Df (Mr))Mr), 

In general, the equation corresponding to S 2m with m > 2 can be written as 
V> m (r + 1) = Df (ip )ip m (T) + Rm(ipo, ih, ■ ■ ■ > VVn-l), (2.24) 

where i? m is a polynomial in (?/>i, . . . , T/> m _i). We stress that the dependence 
on is n °t necessarily polynomial. For example, 

R2 = ^ t 1 D 2 f (^ 1 + f 1 ^ ) 

where ' is used to denote the transposition and f\ is a coefficient of the Taylor 
series f £ = Y. k >o £k fk- 

An analytic solution ■?/>" of the equation f)2.24p can be selected by the 
following asymptotic condition: 

ipm( T ) x i>m{ T ) as r — > oo in the sector (3 < argr < 2n — (3 . (2.25) 

The formal series tp m (r) comes from the re-expansion of the formal separatrix 
given by ( 12. 13ft . In fact, ipm is a formal solution of equation ( 12. 24ft . 

The method of fixing a solution in one region by initial or boundary 
conditions which come from a neighbouring region is known as "complex 
matching ". 

We will now continue with the approximation inside the matching zone, 
we call this domain ^(e): 

T 2 (e) = {r G 7i(e) : Rer > -5 l/ \ -5 1/2 < Imr < 5 1/2 } . 
Theorem 2.6 (Second approximation theorem). 

71-1 

fc=0 

uniformly in T 2 (e). 

We note that this theorem implies the following. 
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Figure 5: Complex matching. 



Matching Zone 




Figure 6: Complex matching. 



17 



Corollary 2.7. There is a parametrisation of the unstable separatrix such 
that 

n-l 

r(r) = J2^k(r) + 0((5rr). (2.26) 

k=0 

uniformly in %.{s). 

Proof. The asymptotic boundary conditions implies V'fc"( T ) = 0(r 2k ). Then 
the theorem with n replaced by 2n implies 

2n-l 
k=0 

71-1 

= ^5 2 V fe -(r) + 0((5r) 2 "). 

k=0 

□ 

2.7 The stable manifold 

So far we have only considered the unstable separatrix. In order to estimate 
the splitting of separatrices we also need to study the stable one. We note 
that the stable manifold is an unstable one for the inverse map. 

Therefore the study of the stable separatrix is analogous to the study of 
the unstable one but with time reversed. The transformation t \— > — t swaps 
the upper and lower half-planes of C and it is convenient to use real-analytic 
symmetry to swap these back. Let 7^ + , T^(e), T 2 {e) C C denote domains 
obtained by reflecting T , Ti(e),T 2 (e) in the imaginary axis respectively. The 
limit flow and approximately interpolating flows are integrable and their 
stable and unstable separatrices coincide. Therefore the same separatrix 
solution arise and the same formal expansion (I2.13P is used for the complex 
matching of the stable separatrix. However, the matching is done on the 
right hand side of the singularity with asymptotic boundary conditions as 

T — > +00. 

Similarly to Theorem 12.41 we conclude that for every n G N there is a 
solution ip + of the equation (II. 5p such that 

= <pf*(t) + 2n+2 , 2n+3 (5) (2.27) 

uniformly for t G (see Figure [7]). Note that we directly started with n 
replaced by 2n. For this solution, instead of (I2.20j) we get 

^+(t) = V f{t) + 2nt2n+1 (r- 1 ) (2.28) 

uniformly in T^~(e) and e G (0, Eq) and instead of Corollary 12.71 we get 
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Figure 7: The shadowed region represents the domain of validity for the local 
approximation theorem for the stable separatrix. 

Corollary 2.8. Uniformly in Zj 7~(e) 

n-l 

^ + (^) = E^ + ( r )+°(( (5r ) 2n ) ( 2 - 29 ) 

fc=0 

where ^( r ) ^ e solution of the asymptotic boundary condition problem 
(12.23p or ( 12.24ft subject to the asymptotic boundary conditions 

V'fe"( r ) x V'&C 7 ") as t ^ oo in the sector —n + /3 < argr < 7r — /3 • 
We note that in general i/jq ^ ip^ [HI [H] . 

2.8 Upper bounds for the splitting of separatrices 

We assume that both tf)~ and ip + are chosen to be closed to (pf 1 . Then 
comparing the estimates ( I2.17P and ( 12.271) 

f(i)-^(f)=0(5 2n ) 

uniformly in T Q n . Then we use the estimates ( I2.28P (with n replaced by 
2n) and (l2~28l) to conclude 

^-(t)-V + W = 2n , 2n+ i(r- 1 ) 

uniformly in Ti(s)r\T^ (e) . Taking into account the difmitions of the domains 
we see that the union of validity domains of these two estimates includes the 
rectangle defined by the inequalities | Im£| < n — and | Rei| < and 
for these values of t 

4) + {t)-4)-{t) = 0{T- n ) = 0{5 n ). (2.30) 
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In the innermost region 7^ (e), closer to the singularity at in, we use a different 
argument. In this region the estimate has to be be close to optimal: we will 
use the upper bound to show that the square of the distance between the 
stable and unstable separatrices is negligible compared to the leading term 
of the splitting (on a line with Imr = — <r log 5). 

The upper bound for the splitting near the singularity is obtained in the 
following way: by Corollaries 12.71 and 12.81 

n-l 

^+(r) - r (r) = £ 5 2k (^ + (r) - V* (r)) + O ((St) 2 ") (2.31) 

fc=0 

uniformly in T 2 (e) PI T 2 + {e). We use that 

ri(T)-^(r) = 0(T m *e- 2 ™ ) (2.32) 

provided —it + [3 < argr < —j3 and Imr < —C2. Consequently 

tp + {r) - Vr(r) = 0(r 4 e" 2 " T ) + 0((5r) 2n ) (2.33) 
provided — c 2 > Imr > — 5~ l l 2 and — ix + (5 < argr < —(3. 

2.9 The flow box 

In order to provide a quantitative description for the difference between tf) + 
and ip~ we use the method based on flow box coordinates. We will prove a 
theorem which essentially says that the restriction of the map f e on a non- 
invariant subset can be analytically interpolated by an autonomous Hamil- 
tonian flow with one degree of freedom. The subset in question contains 
large (but again non-invariant) segments of the stable and unstable separa- 
trices. Simultaneously we introduce "energy-time" coordinates (T, E) for the 
interpolating flow. 

Consider a domain in C 2 defined by 

Us = T x £ = { |Re(T)| < 2, |Im(T)| < - c 3 1 x { \E\ < E (5) } . 

(2.34) 

Theorem 2.9 (The flow box theorem). If c% is sufficiently large, there exists 
a symplectic diffeomorphism S : Us — >• Vs '■= S(Us) such that 

1. f e \ v = where E : Vs — > C is the second component of the map S^ 1 . 

2. S(T, 0) = ip~{T log A<s) (Normalisation). 
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Figure 8: Time component of the flow box coordinates (t := T\og\$). 
3. ||5' _1 || C 2 is bounded uniformly for e G (0,£o)- 

It is convenient to consider S" _1 as a symplectic coordinate map defined on 
Vg- In coordinates (T,E), the Hamiltonian equations with the Hamiltonian 
function E take the form 

T' = 1 , E' = 

and can be easily integrated. The corresponding time-one map has the form 
(T, E) i— > (T + 1, E). It coincides with S^ 1 o $^ o S*. In the original coordi- 
nates the map = f £ . Consequently the first property of S means that S 
conjugates f £ and the translation: 

(T + l,E) = S- 1 of £ oS(T,E). (2.35) 

Let us use (T(x, y), E(x, y)) for the components of S 1-1 and (x(T, E),y(T, E)) 
for the components of S. Equation (12.351) implies 

Tof £ (x,y) = T(x,y) + 1, (2.36) 
Eof E {x,y) = E(x,y). (2.37) 
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By the inverse function theorem, the Jacobian matrix DS 1 = (DS) 1 . In 
coordinate form this equality implies that 



dT 


dT \ 




/ ex 


dx \ 
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dx 


dx 


dy 1 
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dE i 
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0!/ / 




\ dT 


as / 




\ dT 


dT 



where we used that the Jacobian matrix has the unit determinant to get the 
second equality. The second row gives the following Hamiltonian equations: 

dy dE dx dE 

dT = ~~dx~ an df~ = lhj' 

Substituting E = into the equations, using the normalisation condition in 
the form 

(x(T,0),y(T,0)) =^-(T\ogX s ) = (T log X s ) , ^ (T log \ 5 )) 
and switching to derivation with respect to t (where t = TlogA^) we get 

dE 



dx 
dE 
dy 



-log\ 5 ip 2 , 
log X s ipi . 



The last equation has the form of a Hamiltonian system with the Hamiltonian 
E(x, y) after a time rescaling. For future reference we write this identity using 
the vector notation: 

VE(ip-) = - log X 5 #~ (2.38) 
where J = ( J J j is the standard symplectic matrix. 

Remark 2.10. The following facts are not directly used in the proof. It is 
easy to see that equation A2.35\) alone does not define S uniquely. Indeed, 
consider a substitution (T, E) i— > (T, E) of the form 

(T,E) = (T + a(T,E),b(T,E)) 

where a and b are 1-periodic in T. A map S obtained as a result of this 
substitution also satisfy 112. 3 5\) . Although the uniqueness is not important for 
our study, we note that the substitutions of this form exhaustively describe 
the freedom in the definition of S. 

The normalisation property of Theorem \2.9i does not eliminate this free- 
dom completely but only fixes S on the line E = in a way which makes 
comparison between the stable and unstable separatrices more convenient. 
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2.10 Splitting function 

Let us define the splitting function byS 

0(T) = £o^ + (TlogA,5). (2.39) 

We can use to 

1. Find the homoclinic points. 

2. Find the homoclinic invariant uj. 

3. Calculate the lobe areas. 

4. Calculate the distances between ip + and 

Since E — corresponds to the unstable separatrix, zeroes of the splitting 
functions correspond to homoclinic points, i.e., if 0(T"i) = for some 7\ then 
ip + {ti) e W~ for ti = TxlogA. Therefore there is t 2 such that ip + (ti) = 

Let us consider the derivative of the splitting function with respect to T 
at T = T x : 

§ (Ti) = V^l^^miogA) =togA VE%_ it J+( tl ), 

where we used the equality ip + {ti) = ?/>~(t 2 ). We use (I2.38f) to eliminate the 
gradient of E: 

^(n) = - io g 2 A = iog 2 A {ry^ + 

since J* = —J. Taking into account the definition of the symplectic form we 
get 

We remind ourselves that the homoclinic invariant, also known as the Lazutkin 
invariant, is defined by 

Therefore we have established that the homoclinic invariant is proportional 
to the gradient of the splitting function 

u = e'CTi) (2.40) 

log A 



The domain of this function is to be studied. 
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at its zero, 0(?i) = 0. 

Next, we note that O(T) is periodic. Indeed using the definition (12.391) . 
equation (11.51) with t = T log A and the equality (12. 37ft we obtain 

6(T + 1) = £o^+((T + l)logA) = £o/ £ o^+(TlogA) 
= Eo^ + (TlogA) = 0(T). 

Hence we can write the Fourier series 

0(T) = J2®ke 27TikT . 

feez 

Lemma 2.11. Let c\ > be fixed and let p = 7r — c\. If 5 is sufficiently 
small, the domain of the function B contains the strip 

{ T : |ImT| <p/logA 5 }. 

Moreover, for all s G M 

9(s) = Q ie 27ris + eie~ 2wis + 0(e- 4 " p/logAi ) (2.41) 

where 0i = Oi(5) is a first Fourier coefficient ofO. The asymptotic formula 
(I2.4ip can be differentiated with respect to s, the error term does not change 
its order. 

Proof. We note that on the strip the absolute value of 0(T) does not exceed 
Eq(5) defined in (I2.34H . Writing the classical integrals for Fourier coefficients 
of a periodic function and shifting the integration path to the boundary of 
the analyticity strip, we derive 

|©fc| < E (5)e' lklT ^ . (2.42) 

Since log A^ is of the order of S all Fourier coefficients are exponentially small 
(the zero order coefficient being an exception). Moreover, the coefficients 
with \k\ > 2 are negligible compared to the size of the separatrix splitting 
we intend to detect. 

Now we show that Bo is also negligible. We define an auxiliary function 
g(s) = To?/) + (slog A) which describes the T component of the stable separa- 
trix in the flow box coordinates. Using the definition (12.391) . equation (II. 5ft 
with t = T log A and the equality (12.361) we obtain 

g( S + l) = To^ + ((s + l)logA 5 ) =To/ £ o^+( s logA 5 ) 
= To^+(slogA 5 ) + l = g(s) + l. 
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Figure 9: Unstable separatrix is the "x-axis" and the Stable one is a graph 
over it; The area is zero since the map is area-preserving. 

Consequently g' is periodic and its mean value is 1 

/ 1+ (l-g'(s))ds = (2.43) 

for any T\ from its domain. 

Let Ti be a zero of the function 0. Then p 1 = ip + (ti) = ip~{t 2 ) is 
an homoclinic point. The segments of the stable and unstable separatrices 
with their ends at p x and / e (pi) form a closed loop on the plane. The 
total algebraic area of this loop vanishes. Since S is symplectic, this area 
can be evaluated in the flow box coordinates. The segment of the unstable 
separatrix is a straight line which connects the points (T 2 , 0) and (T 2 + 1, 0) 
and therefore the area is given by the integral 

/•Ti+l /-Ti+l 

/ ®{s)g\s)ds = / Q{s)dg(s) = 0. (2.44) 

Now we can come back to estimates for the zero order Fourier coefficient 
of the function which is given by the integral 

O = / Q(s)ds. 

Taking into account (12.441) we rewrite it 

O = [ Tl+1 Q(s)(l-g'(s))ds. 

JTi 
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Then using equation (I2.43P and the fact that O is constant we get: 



rTx+l 



6 



o = 



/ (e(s)-e )(i- g >(s))ds 



Jti 



Under the integral we see the product of two functions, each one is periodic 
and of zero mean. The Fourier series arguments show that 



e(s) - 6 = 0(e- 2np/logXs ) and 1 - g'(s) = 0(e~ 2wp/l ° sXs ) (2.45) 



as the exponent in the right hand side is almost double the expected leading 
order of the separatrices splitting. 

In particular, combining the last upper bound with (12.421) we conclude 
that on the real axis 



Moreover, this function is almost sinusoidal since if we keep apart the two 
largest Fourier coefficients the sum of all others is very small and for s G M 
we get (EH}. 

Of course the Fourier coefficients 0i and 6_i = depend on 5 and are 
exponentially small due to the upper bound (I2.42p with \k\ = 1. Nevertheless 
the the error term is much smaller than the upper bound for 0±i. In the 
next section we will construct an asymptotic expansion for B_i and conclude 
that generically it really dominates the error. □ 

Finally we show how to use the splitting function to estimate the lobe 
area between the stable and unstable manifolds. Suppose 7\ and T 2 are two 
zeroes of G. Then the lobe area is represented by the following integral: 



Taking into account the upper bound (I2.45P and (I2.46P for the functions 
under the second integral we conclude that the lobe area equals to 



since s is real. Consequently, Qq is extremely small: 



e = 0(e- 4np/logXs ), 




(2.46) 
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Figure 10: Domain used to estimate the first Fourier coefficient of the split- 
ting function 

2.11 Asymptotic expansion for the homoclinic invari- 
ant 

Consider a line in the complex plane of the variable T defined by 
£(5) = {ImT = Pl (6) , |ReT|<l}, 

where 

log As 

and a = - is a positive constant. We fix the following relation between the 
parameters: t = T log \$ and 

t — in ^ in . d . 

log X s log X s 

We estimate the splitting function defined in (12.391) by expanding the energy 
E in Taylor series centred at a point ip~ of the unstable manifold. Using E o 
ip~ = 0, equation (I2.38I) for the gradient of E and the uniform boundedness 
of the second derivatives of E we find that 

0(T) = logA 5 det (^t--^--^\ +0(|^--^ + | 2 ) . 

The functions in the right hand side are evaluated at a point r with Im r = 
—a log . Using (12.311 with n replaced by 2n gives us the following estimate 
for the difference of ip^ 

2n-l 

^ + - r = E 62k ^t - + oc^iog^r 1 ). 

k=0 



27 



We note that on £(5) 

by the upper bound (12.321) . hence it follows that 



(2.48) 



n-l 



lf) + — if) 
lf) + — if)' 



S 2 \^t - ^k ) + 0(5 4n log" 1 (T 1 ) and 



k=0 



0(S 2n log m 5~ u 



dr 



< 1 



where m = max{4n, mo, . . . , m2„_i}. On a sufficiently large neighbourhood 
of £(S) we have |^~(r)| < 1, therefore the Cauchy estimate implies 
and consequently 

'dip' 



0(T) 



det 



dr 

dtp" 



ip~ -ip + ) +0{\ip- -ip + \ 2 ) 



n-l 



^ I ^pX)**^* ] +0(5 4 "log 2m 5- ] < 

fc=0 



Now, since if) is analytic the series (I2.26j) can be differentiated term-wise, 
so using (12.481) we obtain 



'n-l 



6(T) =det J2 S 



,fc=0 



dr 



0(5 4n log" 1 ^ 1 ), 



k=0 



where m = max {2m, (2n — l)m } . We use that the determinant is bi-linear 
to rewrite this estimate in the following form 



n-l ri— 1 



fc 1= 0fc 2 =0 



( T ) = EE 52(kl+k2) det ( -JT> K - < ) + 0(5 An log™ r 1 ). (2.49) 



Next we estimate the derivative of if>^ . Since it is an analytic function in 
a small domain T> we can use the Cauchy estimate to estimate it using if>^ 
itself. Using that if>^ is asymptotic to, see (12.251) . the series given by (I2.12p 
yields 

d ^ kl - 0(r mfc i) = 0(log mfc i <T lN 



dr 



which together with the estimate (12.481) gives 



dr 
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The last estimate shows that terms in (12.491) with k% + k 2 > n are not larger 
than the error term. We collect terms with {k\ + k 2 ) < n to get 

6(T) = |V J2 det^^^-^+O^log^r 1 ), (2.50) 



k=0 ki+k2=k 
ki,k 2 >0 



where fh = max kl+k2=n {m,m kl + m k , 2 }. 
Let us introduce the auxiliary functions 



6 k (r)= det(^(r);^ 2 (r)-^ 2 (r)y (2.51) 

k 1 +k 2 =k ^ T ' 

We note that these functions are independent from the parameter 5. The 
next lemma gives a more accurate estimate for these functions. 

Lemma 2.12. For every integer k > there exist u k G C and iVfc G N such 
that 

6 k (r) = u k e- 27riT + 0(r N *e- A7riT ), (2.52) 
where r G {— 7r + < argr < — /3 1; Imr < c}. 
According to the lemma we have 

fc ( r ) = c^e" 2 ™" + 0{5 An \og Nk cT 1 ) (2.53) 
on the line £(S). Substituting this and (12.511) into equation (I2.50p . we obtain 

n-l 

6(T) = S 2k u k e- 2mT + 0{5 4n log M r 1 ), 

k=0 

where M = maxjm, iVo}. Now we can evaluate the first Fourier coefficients 
of 6: 

rl+ipi(8) 

9_i = / e 2nlT Q(T)dT 

Ji Pl (8) 

/1-ialogS- 1 /n-l \ 

e 2 ™ ^W" 2 ™ + O (5 4n log M (T 1 ) dr , 
-«o-1ob;<5 _1 \ ,__„ / 



e log a 4 



vifc=0 



where the exponential factor comes from the change of the variables ( 12.471) . 
Expanding the parenthesis and taking into account that \e 2mT \ = 5~ 2n on the 
integration path we get 



'n-l 



log A 



,fc=0 
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With this bound on the first Fourier coefficient we are now ready to estimate 
the size of 0(i) for t real. In particular, if s G K then (I2.4ip with p = |7r 
implies 



9(s) = Q^e~ lms + G^e 2 ™ + O [ 
n-l 



2tt2 _ / _ _ 27T 2 



fe=0 ^ 

There are two cases to consider. If Uk vanish for all k, this formula simply 
gives an upper bound for the derivative of 9 and, as a result, for the homo- 
clinic invariant of any primary homoclinic trajectory. If, on the other hand, 
some 00^ ^ the leading term is larger than the error. In this case it is more 
convenient to rewrite the formula in the form 

n_1 / n_1 \ -2 2 / 2 \ 

0(s) = a k$ 2k cos 2vrs + ^S 2k e^r + q I &» e -&x J , 

fc=0 V fc=0 / ^ ' 

where a® = \uq\ and at least one of the amplitudes does not vanish. The 
implicit function theorem then implies that the function has exactly two 

zeroes per period. The derivative of O at the zeroes is ±27re lo s A X]fc=o a k$ 2k - 
Therefore there are exactly two primary homoclinic orbits and the relation 
(12.40p implies our main result: 



27T " 



CO 



"log A 2 

ta fc=0 



a k 5 2k e~^ + O I 5 2n e~^ J 

h— n \ / 



3 Normal form for the bifurcation 

3.1 Formal series and quasi-homogeneous polynomials 

In this section we will mainly be interested in transformations given in the 
form of formal power series. We will consider the series in powers of the 
space variables (x, y) combined with expansions in the parameter e. The 
series have the form 

9{x,y,e) = c ki m x k y l e m . 

k,l,m 

A formal series is treated as a collection of coefficients. Formal series form an 
infinite dimensional vector space. Addition, multiplication, integration and 
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differentiation are defined in a way compatible with the common definition 
on the subset of convergent series. 

The series involve several variables and it is convenient to group terms 
which are "of the same order". A usual choice is to consider x and y to be 
of the same order. But for purpose of this paper it is much more convenient 
to assume 

• x is of order 2; 

• y is of order 3; 

• e is of order 4. 

Then a monomial x k y l e m is considered to be of order 2k + 3/ + 4m. We can 
write 

g(x,y,e) = ^2g p (x,y,e) where g p (x,y,e)= ^ c Mm x k y l e m 

p>0 2k+3l+4m=p 

is a quasi-homogeneous polynomial of order p. We stress that this notation 
does not refer to a resummation of the divergent series but simply indicates 
the order in which the coefficients are to be treated. 

In order to give a rigorous background for manipulation with formal series 
we define a metric on the space of formal series Sj. Let g and g be two formal 
series and let p denote the lowest (quasi-homogeneous) order of g — g. If 
g 7^ g then p is finite and we let 

d(g,~g) = 2- p , 

otherwise we assume 

d(g,g)=0. 

It is a straightforward to check that (S), d) is a complete metric space. More- 
over polynomials are dense in (fj, d). Hence we can define formal convergence 
and formal continuity on the space of formal series. In particular an operator 
is formally continuous if each coefficient of a series in its image is a function 
of a finite number of coefficients of a series in its argument. 

Let x an d g be two formal power series. We note that any of the series 
involved in next definitions may diverge. The linear operator defined by the 
formula 

L x9 = {g,x} (3-1) 

is called the Lie derivative generated by x- We note that if x starts from 
an order p and g starts with an order q, then the series L x g starts with the 
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order p + q — 5 as the Poisson bracket involves differentiation with respect to 
x and y. If p > 6 the lowest order in L x g is at least q + 1. 
We define the exponent of L x by 

exp(L x )s = ^L^, (3.2) 

fc>0 

where L x stands for the operator L x applied k times. The lowest order in 
the series L x g is at least q + k and consequently every coefficient of exp(L x )g 
depends only on a finite number of coefficients of \ and g. 

If the lowest order of x is a t least 6 the series (13.21) converges with respect 
to the metric d, i.e., each coefficient of the result is function of finite number 
of coefficients of x an d <?• 

We consider consider the formal series 

$ x = (exp(L x )x,exp(L x )?/) . 

We say that $ x is a Lie series generated by the formal Hamiltonian x- If 
X is a polynomial the series converge on a poly-disk and coincide with a 
map which shifts points along trajectories of the Hamiltonian system with 
Hamiltonian function \- F° r this reason sometimes we will call this series a 
time one map of the Hamiltonian system x- 

We note that it is easy to construct the formal series for the inverse map: 

^x 1 = (exp(-L x )a;, exp(-L x )y) . 

Then <3> x o $ x 1 (x,y) = (x, y). We also note that 

go$ x = exp(L x )g. 

These formulae are well known to be valid for convergent series and can be 
extended onto S) due to the density property. 



3.2 Formal interpolation 

We study the following family of area-preserving maps 



x + y + f(x,y,e) 

y + g(x,y,e) 



and the Taylor series of f(x, y, 0) and g(x, y, 0) do not contain constant and 
linear terms. Since F e is area-preserving we have 



det DF F 



1+dJ 1+dyf 

d x g 1 + d y g 
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which is equivalent to 

d x f + d y g + {f,g}-d x g = 0. (3.3) 

Next theorem states that F £ can be formally interpolating by an autonomous 
Hamiltonian flow. 

Theorem 3.1. Let F £ be family of area-preserving maps such that 

DF (0)=(J \ 

then there exists a unique (up to adding a formal series in powers e only) 
formal Hamiltonian h e such that 

*JL = F <- ( 3 - 4 ) 

Proof. The Taylor series for F £ : (x,y) i— > is written as a sum of 

quasi-homogeneous polynomials: 



x x =x + y+ f P (x,y,e), 

p>4 

yi = y+ J2g P ( x ,y, £ ), 

p>4 

where f p and g p are quasi-homogeneous polynomials: 



(3.5) 



f p {x,y,e)= £ fki m x k y l e m , 

2k+3l+4m=p 

g p (x,y,e)= £ gki m x k y l e m . 

2k+3l+4m=p 

Similarly we look for h £ in the form of a sum of quasi- homogeneous polyno- 
mials using the same quasi-homogeneous ordering for terms: 

h £ (x,y) = ^h p (x,y,e) where h p (x,y,e)= h Hm x k y l e m . 

p>6 2k+3l+4m=p 

The time-1 map of this hamiltonian is given by the Lie series 



x + L he x + h L hX 

k>2 

1 rk 



y + L he y+j:±L k he y 

k>2 



where the operator L he is defined by 

dtp dh £ dip dh £ 



L h e (tfi) = W,h £ } 



dx dy dy dx 
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Now we use induction to show that there is a formal Hamiltonian h e such that 
= F £ at all orders. This is an equality between two formal series, i.e., all 
coefficients of two series coincide. We note the first component of the series 
starts with the quasi-homogeneous order 2 and the second one starts with 
3. Therefore it is convenient to consider an order p in the first component 
simultaneously with the order p + 1 in the second one. In this situation we 
say that we consider a term of the order (p,p+ 1). 

Equality (13.41) is equivalent to the following infinite system 



k>2 ' k 

K+a2/ + E^[ L t4 = 9p+i, (3-7) 

k>2 

where [-] p is used to denote terms of the quasi-homogeneous order p of a 
formal series. We note that 

order(L/ lfc x) = k — 3 and order(L hk y) = k — 2 . 

First we check that the equations can be solved for p = 3. It is easy to see 
that 

dh 6 



L h 6 y 



dy ' 
dh 6 



dx 

Then the equations (I3.6p . (13.71) take the form 



y, (3i 



dy 
dh H 



dx 



Since is a quasihomogeneous polynomial of order 4 it is independent from 
y which is of order 3. Consequently 

g 4 (x,y,e) = g 4 (x,e) . 

Then system (13.81) has a unique solution of the form 

y 2 

h = y + G 6 (x,e), 



34 



where G 6 is a quasi- homogeneous polynomial of order 6 in re and e. 

Let us now proceed with the induction step. Assume we have h 6 , . . . , h p+2 
and want to compute h p+3 for p > 3. To make the combinatorics easier we 
introduce the notation 

L s ip = {if, h s+5 }. 

Then 

Lh e ¥ = W, M = fa> M = ^2 L s¥- 

p>6 s>l 

Note that L s maps a quasi-homogeneous polynomial of order j into a quasi- 
homogeneous polynomial of order j + s. So we have 

k 



^3 = \5l Ls j W 

= E E 

m>k si+S2+-..+s k =m 



E 



si+s 2 +...+s k >k 

Sl,S2,...,Sfc>l 

L S1 -- ■ L Sk ipj. 



Now let us consider the components of at order (p,p + 1), the first 
component is 



L ^ x + E ^ 



fc>2 



X 



P-2 x 

l ^ + Ett E 



and the second component 



fc>2 



fc! 

k=2 si+S2+...s k =p-2 



p-2 



L heV + J2^ L iy =L p- 2 ^ + E^! E L Sl ...L Sk y. 



p+1 



We want 



We have 



2 nd comp 



k\ 

k=2 s\+S2+---s k =p-2 



u 

9p+i 



{x, h p+3 } , 

{y, h P +?>} ■ 



Thus 



/ dh p+3 \ 








1 dhp+3 i 









p-2 ^ 
/p ~~ E fc! 



E 

fc=2 si + ...+Sfc=p-2 
p-2 



L S1 . . . L Sk x 



\ 



, 9p+i - E F E L S1 ... L Sk y 

\ k=2 s 1 +...+s k =p-2 / 



(3.9) 
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The Hamiltonian h p+s is defined from this equation uniquely up to a function 
of £ if and only if the right hand side is divergence free. The last condition 
follows from the area-preservation property due to the following lemma. 

Lemma 3.2. If expansions of two area-preserving maps coincide up to the 
order (p — l,p) then the divergences of the terms of the order (p,p + 1) are 
equal. 

Proof. We use the area-preserving property ( 13.31) and collect the terms of 
order p — 2 

d x f p + d y g p+1 - d x g P + ^ {/*, gi} = 0. 

k+l=p+3 
k,l>4 

Consequently, 

d x fp + 9 y g p+ i = d x g p 

= d x g p 

□ 

We finish the proof by the following observation. Let h = X]fc=6 -^ or 
any choice of h p+ s the maps $| and F £ satisfy the assumptions of the previous 
lemma therefore their orders (p,p+ 1) have the same divergence. That is the 
solvability condition for the equation (I3.9p . □ 

3.3 Simplification of the interpolating Hamiltonian 

The Hamiltonian h £ can be simplified with the help of a formal canonical 
change of variables. 

Theorem 3.3. If h(x,y,e) = J2 p >6^p( x ^y^ £ ) such that the leading order 
has the form 

he(x,y,e) = — + ay - hex (3.10) 

then there exists a formal canonical substitution such that the Hamiltonian 
takes the form 

~ y 2 

h(x, y, e) = — + ^2 u p( x ' £ )- 

p>6 



k+l=p+3 
fej>4 



Y {fk,9i\ = 9 x f P + d y g p - 



k+l=p+3 
fcj>4 
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Proof. We prove the theorem by induction. The order 6 is already in the 
desired form, just let w 6 = ax 3 / 3 — hex. Then suppose we transformed the 
Hamiltonian to the desired form up to the order p with p > 6. For any 
quasi-homogeneous polynomial Xp °f order p we have 

h = e L *p h = h + {h, x P } + O 2 p_i +6 
= h+ {h 6 , x P } + P +2 

We see that this change does not affect terms of order less or equal p. Col- 
lecting all terms of order p + 1 we get 

hp+i = h p+1 + {h,x P } ■ (3.11) 

In order to complete the proof we show that there is a polynomial Xp such 
that h p+ i is independent from y. The situation is different for p even and p 
odd. We will show that if p is even then for any h p+ i there is Xp such that 

h p+ i + {h 6 , Xp } = 0. (3.12) 

The polynomial Xp is defined up to addition of an arbitrary polynomial of he 
and e. On the other hand, if p is odd then there are two unique polynomials 

u p+ i = u p+1 (x,e) and x P +i = Xp+i(x,y,e) such that 

hp+i + {h 6 , Xp} = u p+ i . (3.13) 

In order to prove these statements we note that since y is counted as a 
term of order three and terms of order two and four respectively, a 

quasihomogeneous polynomial of even order contains only even powers of y 
and one of an odd order contains only odd powers. Let p + 1 = 3q + r with 
r G {0, 1, 2}. The largest possible power of y in h p+1 does not exceed q and 
we can write 

q/2 (5-l)/2 

h P+i = ^2y q ' 2ls r+ei(x,e) and Xp = ^2 y q ^ 1(J r+&i+2{.x,e), 

1=0 1=0 

where Sj and Oj are quasi-homogeneous polynomials of x and e. Taking into 
account 

du 6 dxp dx P 

{h,Xp} = — o a ?/^— 

ax ay ox 

and collecting terms in h p+ i + {/ic Xp} which have of the same order in y 
(excluding the y-independent terms) we get 

-a— h s r = , 

-a — ^ h (g - 2/ + 1)— a r+6 i_ 4 + s r +6« = 0. 
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The first equation defines cr r+2 . If r ^ 2 the solution is unique, otherwise a 
constant times e can be added. Then oy +w+2 are to be chosen recursively for 
1 < / < g/2 to satisfy the second equation. 

The theorem follows by induction in p. □ 



4 Formal separatrix 

In this section we will describe properties of an asymptotic expansion for the 
separatrices providing a proof for the statements of Section 12.21 



4.1 Auxiliary functions 

In this subsection we describe a useful class of functions. A function belongs 
to this class if 

1. it is periodic with the imaginary period 2ixi. 

2. it is analytic on the entire C except for poles at 7ti(2k + 1), k G Z. 

3. it vanishes as Ret —>■ ±oo. 

Any function which satisfies these three assumptions is a polynomial, without 
a constant term, of two "base" functional 

* = ^I m = -^I- (41) 

Equivalently, the function can be written in the form P(t)q) +i]iQ(i]q) where 
P, Q are polynomials of 7] only and P(0) = 0. Indeed, a straightforward 
substitution shows that r/ , r/i satisfy the differential equation 

Vo = Vi, (4-2) 
3 2 

Vi = Vo-^Vo- 
This equation is Hamiltonian with the Hamiltonian function 



2 We note that tanh | is 2m periodic and analytic on C except for simple poles at 
i7r(2fc + 1) with k € Z. Suppose a function (p(t) satisfies properties 1 and 2 and is bounded 
as Re t — > ±oo. Comparing Laurent expansions around any of the poles, we can construct a 
polynomial p such that tp(t)—p(tanh |) has no singularities and, consequently, is constant. 
The condition 3 adds a restriction on coefficients of the polynomial (p(l) = p(— 1) = 0). 
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The Hamiltonian is constant along solutions of the differential equation. 
Since both 770 and r]i vanish as Re — > 00, we get the identity 

v!-ri 2 +Vo = Q- (4-3) 
This last equality and (14.21) imply that 

V 2 o = Vl = V 2 o ' Vo • (4-4) 

Using this identity we can exclude all even powers of 771 and reduce all odd 
powers to the power one. These identities will be used to simplify formulae 
involved in the construction of the formal separatrix. 



4.2 Formal separatrix of the flow 

In this section we prove Theorem 12. 1[ We remind that S = e 1 ^ 4 . By Theorem 
13.31 there is a formal canonical change of variables which transforms the 
Hamiltonian h £ to the simpler form 

H e (x,y) = ^ + U(x,e), (4.5) 

where 

U(x,e) = J2 u kmx k e m . (4.6) 

m>0,fc>l 
fc+2m>3 

After this change equation (12.71) . which we for convenience repeat 

^(<J)X= JV# £ (X), (4.7) 

takes the form 



/i(e)x = y, 
H{e)y = -U'(x,e) 



(4.t 



where U' = This system can be conveniently rewritten in an equivalent 
form as a scalar equation of second order: 

(fi(e)) 2 x = -U'(x,e). 

In order to simplify notation we introduce the new auxiliary series 

Since \i was an odd series in S, the series (5 is even. The equation takes the 
form 

2(3(e)x = -U'(x,e) 
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Multiplying with x and integrating with respect to t we obtain 

(3(e)x 2 = -U(x,e) +c(e), (4.9) 

where c(e) is a constant of integration. Now we use the following ansatz for 
solving this equation: 



X(t,£) = ^5 2fc X fc (t), 
fe>l 

(3(e) = ]Ta fc 5 2 \ 



fe>i 



E c ^ 2fc - 



and 

c(e) 

k>3 

Substituting these series into equation ( 14.91) we get 

E^ 2fc ) (e^w) + E «^fe*W) ^ 4m = E c ^ 2fc - 

\fc>l / \fc>l / k+2m>3 \j>l / fc>3 

(4-10) 

We recall that this equation is considered in the class of formal series which 
assumes that two formal series are equal if and only if their coefficients co- 
incide. We will use induction to show that the equation can be satisfied at 
every order in 5 2 . 

We note that in (I4.10p the least order of 5 is 6. Collecting all terms of 
this order we get the following equation 

a x x\ + u m x\ + U\\x x = c 3 . (4.11) 

This equation looks very similar to (14. 4p but has different coefficients in front 
of its terms. This suggests x\ to be of the form 

x 1 = b + b 1 r ]0 , (4.12) 

where bo, b\ £R are to be determined. We insert x\ into (14. lift which gives 

aiihrjo) 2 + u 30 (b + b^f + u n (b + b^o) = c 3 . 

Using equation (14.41) to replace ?) 2 and expanding the parentheses we rewrite 
this expression as follows 

c 3 = a^l (r)l - 7]l) + u 30 bl + u m d,blbi7]Q + w 30 3&o& 2 ?7o ( 4 -l 3 ) 
+u^b\rjl + u n b + UnbtfjQ. 
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We now collect terms in (14.131) of equal powers in r] and obtain the following 
system of equations 

-a x b\ + u 30 bl = 0, 



axb\ + u 30 3b bl = 0, 



(4.14) 



k c 3 = u u b + u 30 bl. 
The last equation defines c 3 . The first three lines imply 



b 



± 



3«30 : 



-36o, 



(4.15) 



ai = -3u 30 b . 



We make the following two remarks. Note that, to the leading order, bo is 
the x coordinate of the equilibrium point and c 3 is its energy. 

The sign of &o should be chosen to ensure a± > since the latter coefficient 
is the first term of the expansion of /3(e) which is a square of a real series. 

We have seen that equation (I4.1UI) has a solution at the leading order <5 6 , 
and continue by induction. We explain the first step of the induction in more 
detail in order to illustrate the method used in the general step. 

Let us collect the terms of order 5 8 which gives the equation 

2a 1 x 1 x 2 + a 2 x\ + u n x 2 + u 2 \x\ + 3u 30 xjx 2 + u 40 x^ = c 4 . (4.16) 

We note that this is a linear non-homogeneous equation in x 2 . Therefore it 
defines x 2 uniquely up to addition of a solution of the corresponding homo- 
geneous equation. 

Using the form of X\ given in equation f)4.12p and the relationship (14. 4|) 
we rearrange this equation as 



c 4 = 2(XiXiX 2 + (3m 3 o2i + Uu) x 2 

+a 2 bj(r]l - r\l) + m 2 i(& + Mo) 2 + M 4o(&o 



(4.17) 



Mo) z 



Before solving this equation for each power of 770 separately we make some 
observations which will simplify our analysis. First we note that the second 
line of (14.171) is a polynomial in r) Q of order 4. The first line of (I4.17P tells 
us that x 2 should have a pole of order 4, hence we assume that x 2 takes the 
following form 

b 20 + b 21 r] + b 22 r]Q. 



x 2 
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Now let us simplify the first line of (14.171) beginning with the second term 

(3u 30 xl + Wii) = 3u 30 (b + Mo) 2 + «n 

= 3u 30 bl + 6ti 3 o&oMo + 3m3o&i^o + u n- 

Using (j4.15p we conclude that 

(3u 30 xl + = 3u 30 bib (2r] - 3^) (4.18) 
To ease the notation we let 

A = 3u 3 o&i&o, (4.19) 

and we notice that A 7^ by virtue of ( 14.15p . Next we simplify the first term 
in the first line of (I4.17P 

Xx± 2 = (621 + 2622^0) VohVo 

= 61 (621 + 2622^0) (Vo ~ Vo) , 

where in the last step we used (I4.4p . Then equation (14.171) reads 

c 4 = 2axh (621 + 2b 22Vo ) {vl ~ Vo) (4-20) 
+A(2r] - 3?7o) (b 20 + 621^0 + b 22 rjl) 
+a 2 b\(r]l - t]q) + u 21 (b + b^) 2 + w 4 o(6 + Mo) 4 - 

Let us now solve (I4.20p for x 2 ,c^ and a 2 . The strategy is to solve (I4.20p 
in each power of 770 separately in the following order 

1. C4 is defined from the order in 770. 

2. 620 is defined from the order 1 in rjo. 

3. 622 is defined from the order 4 in 770- 

4. 621 and a 2 solve a system of equations obtained from the orders 2 and 
3 in rj . 

At order in r] we get 

c 4 = u 21 b 2 Q + u 40 i>o, ( 4 -21) 

which defines the coefficient C4. Next we proceed with the terms of order 1 
in ?7o 

2Ab 2Q + m 2 i26(A + «4o4&o&i = 0, 
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which gives 

u 2 i2b bi + u 4 o4&o&i . . 

°2o — ^4 • I 4 - 22 ) 

Since A is nonzero equation (I4.22p defines the coefficient 620 uniquely. At 
order 4 in 1] we have 



hence 



-(4tti&i + 3A)&22 + = °> 



^22 



4ai6i + 

The denominator in (14. 23ft is non-zero since, using (14. 15ft and (14.191) , 



(4.23) 



Aa\bi + 3A = -12u 3 ob bi + 9m 30 6o&i = 9ii 30 &f 



0: 



and U30 7^ by assumption and 60 7^ by ( 14.15B . Therefore (I4.23j) determines 
622 uniquely. The terms of order 3 in rj are 

-a 2 b\ + Au m b b\ + (-2ai&i - 3A)6 2 i + ^a x b x b 22 + 2A6 22 = 0, 

and the terms of order 2 in r] are 

-3A6 20 + 2a 1 b 1 b 21 + 2A6 2 i + a 2 &i + u 2 ibj + W406&0&? = 0. 

These two terms define the following system of equations in 6 2 i an d «2 

2a 1 fc 1 + 2A b\ \f b 21 \ = ( 3Ab 20 - u 21 b\ - u 40 6&§&? 
-2ai&i - 3A -6f y y a 2 / \ -4a 1 6 1 6 22 - 2A6 22 - u m Ab b\ 

(4.24) 

Call the matrix in the left-hand side M. This system of equations has a 
unique solution since 

det M = Ab\ = -8lu 30 b* ^ 0. 

Therefore 621 and a 2 is the unique solution of (I4.24p . We have now shown 
that x 2 ,a 2 and C4 are uniquely defined. 

Now let us continue with the inductive step at an arbitrary order of 5 2 . 
For a general n the terms corresponding to 5 2n are 



k+l+m=n k+2m>3 ji+...+jk=n—2m 

h,l,m>\ ii>-,ifc>l 



•^31 ■ ■ ■ ^Jk ~ Ui ri 



This equation has the form 

2aiXi± n + (3u m x\ + u n )x n + P{xi, . . . , x n -i) + a n x\ = c n , (4.25) 
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where P(x±, . . . , x n _i) is a polynomial. We note that P is a (known) polyno- 
mial of order n+2 in tjq and that all coefficients are unique since our induction 
assumption is that x±, . . . ,x n -i are unique. We will show by induction that 
there exist a unique a n and a unique 



» = bnmT lo- (4.26) 



m=0 

That is, we will show that a n and all b nm are uniquely defined. 
Using equality (14. 181) which states that 

(3u 30 xl + Mn) = A(2t] - 3?7o), 

where A = 3W30&1&0 7^ (see equation (14.191) ) and that x n has the form (14.261) 
we simplify (14.251) as follows 

c n = 2a x b x {r]l - rfc) (b nl + 2b n2 r] + ... nb^- 1 ) (4.27) 
+A(2t 1o - 3 V l) (b n0 + ... + b nn ri%) 
+a n b\{ril - r)l) + P{x x , x„_i). 

We will now solve this equation for each power in t]q separately. The strategy 
is the following 

1. The power will define c n uniquely. 

2. The power 1 will define b n Q uniquely. 

3. The power n + 2 will define b nn uniquely. 

4. The powers j G {4, . . . , n + 1} (starting from j — n+ 1 and proceeding 
in decreasing order) will define the coefficients & nj -2 uniquely. 

5. Finally, at powers 2 and 3 we will obtain a system of equations in b n i 
and a n which we show has a unique solution. 



In what follows we use 



P 



to denote the terms in the polynomial P(x\, . . . , x n -i) 

k 

which are of order k in t]q. Let us now start solving (14.271) at order in r/o- 
At this order we have 

c = P 

L J 

which defines c n uniquely. Next, at order 1 in t]q we have 

2Ab n0 + \P] = 0, 
L J l 
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thus 



b n .O 



p 

- Jj 

2A 



Since A ^ this equation defines b n0 uniquely. 
At order n + 2 in i] Q we get 



thus 



-2a\nb\b r , 



b„ 



3Ab nn + 



n+2 



P 



n+2 



2na 1 b 1 + 3A 

By equation (14.151) g^Oi = A, thus the denominator is equal to A(2n + 3) 
which is non-zero, hence b nn is defined uniquely. 

Now we continue with the orders j £ {4, . . . , n + 1}. At order j we solve 
the equation 



(2oi6i0' - 2) + 3A) - 6 rw -_i(2oi6i(j - 1) + 2 A) - 



P 



0. 



thus 



6 nj -_i(2ai6 1 (j-l) + 2A) + 



P 



2ai6i(j -2) + 3A 

Since we are proceeding in decreasing j the coefficient in the right-hand 

side has already been determined. The denominator 

2a A (j - 2) + 3 A = A(2(j - 2) + 3) 

is non-zero since j — 2 is positive for all j £ {4, . . . , n + 1} . Therefore the 
coefficients & n j-2 are uniquely determined for j £ {4, . . . , n + 1}. 

Finally we have to solve for the remaining coefficients a n and b n ±. They 
are obtained from the equations at order 2 and 3 which are given by 



2ai6i&„i + 2Ab ni - 3Ab n0 + a n b\ 



P 







and 



4ai6ife n2 - 2a 1 b 1 b nl + 2Ab n2 - 3Ab nl - a n b\ + 



P 



0. 



respectively. We observe that this system of equations may be written as 



2a x 6i + 2A b\ 
-2ai6i - 3A -b\ 



bnl 

a„ 



3Ab, 



nO 



P 



J 2 



-4ai6i6 n2 - 2A6 



'n2 



P 



J 3 



(4.28) 
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The determinant of the matrix in the left-hand side is Ab\ = — 81ii3o&o which 
is non-zero, hence b n i and a n is the unique solution of (14.281) . 

We have now shown that 6 n0 , . . . , b nn and a n and c n are uniquely defined, 
hence x n is a unique solution of (14.25!) . 

It follows by induction that x m is uniquely defined for all m. 

Now we can reconstruct the second component of the solution of sys- 
tem (14. 8p . Since 

x = J2$ 2k Mt) 

k>l 

the second component of the solution is restored by 



k>l 



k>l 



Since Xk(t) is a polynomial in 7]o of order k we conclude that the solution of 
the system H4.8[) can be written in the form 



( x= Y J 5 2k x k {t), 



k>l 



(4.29) 



k>l 



where yu-i is a polynomial in t]q of order k — 1. 

The solution (I4.29[) is a formal solution of the system of equations (14.81) . 
Equation (14 .8p was obtained from (12. 7p by a change of variables 



i(X) 



(4.30) 



where Xe=E p >6Xp where Xp is a quasi-homogeneous polynomial of order p. 
To finish the proof of the theorem we have to invert this change of variables 
and consider the structure of the resulting formal series X. As the change 
of variables is close-to-identity (in the quasi-homogeneous sense) it is easy to 
invert (I4.30p . For our purposes it is sufficient to consider the change in its 
most general form, hence 



X 



/ x+ £ d klm x k y l 6 4m \ 

2k+3l+4m>3 

, y+ E fki m x k y l 5 4m , 

\ 2fc+3i+4m>4 / 



(4.3i; 



where dki m , fkim are coefficients and (x,y) are given by (14.291) 
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Let us now consider the structure of X, that is, let us consider the sums 
in the right hand side of the previous equation. Obviously 5 Am is an even 
power of 5. We also note that any power of the series x is a series of the same 
type as x itself, i.e. a series in even powers of 5. This is not the case with y . 
An even power of the series y gives a series in even powers of 6 while an odd 
power of y gives a series in odd powers of 5. Hence x k y l 5 4m is a series in even 
powers of 5 if / is even and in odd powers if / is odd. 

Furthermore, since (x,y) are series of polynomials in rjo(t) and r)o(t) any 
power of them are series of polynomials as well. Hence we conclude that X 
is a power series in both odd and even powers of 5 and that each power of 5 
is multiplied by a polynomial in r] (t) and r]o(t). 

A more careful study reveals even more about the structure. Consider 
( I4.29P again. We see that taking any power of x gives a series of terms of the 
type 

5 2k P k (t), 

where is a quasi-homogeneous polynomial in r] of order k. When taking 
powers of y we note that 



.fe _ • k 



(Vo ~ ?7o) fc//2 > V even 



Vl ~ Vo ~ I ?7o(% 2 - r/g)^- 1 )/ 2 , p odd • 
Bearing this in mind an even power of y gives a series with terms of the type 

S 2h Q<k(t), 

where Q<k{t) is a homogeneous polynomial in ?7o(t) of order at most k. Sim- 
ilarly, taking an odd power of y gives a series with terms of the type 

5 2k+1 Q< k -i(t), 

where Q<k-i{t) is a polynomial in ?7o(t) of order at most k — 1. Based on 
these observations we conclude that the form of the formal separatrix X is 



X(t,e) 



( E^4 + ^WE5 2p+1 4-i \ 

\ P>1 P>2 / 



(4.32) 



where x^x^y^y^ denotes polynomials in ^o(^) °f order p. This finishes the 
proof of the theorem. 
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4.3 Re-expansion near the singularity 

In this section we will derive equation ( 12. lip , i.e. we will expand the formal 
separatrix X in its Laurent series around the singularity t = in. We do this 
by first expanding the two "base" functions r)o(t) and T]o(t) into their Laurent 
series around the singularity. The Laurent series of r}o(t) a t t = in is given 
by 



r] {t - in) = ~ ( 1 + E c kt 



2k 

k=l 



where Ck G M V£j G N. The distance to the next singularity is 2n, which is 
then the radius of convergence for this series. Secondly we note that the 
Laurent expansion of r)i (i) at t = in is given by 



r/i (t - iir) = — ( 1 + d k t 



2k 

fc=l 



where d k G M VA; G N and it converges for \t\ < 2n. 

We are now ready to determine the Laurent expansion of (14.321) . Let us 
consider its first component which we write as 



n-l 



[x]„-i = E^w ( 4 - 33 ) 



k=2 



where ip\{t) is the term corresponding to S h in the formal separatrix ( 14.321) . 
ip\ has a pole of order k. Using the Laurent expansions of T]o(t) and fji{t) 
respectively, the Laurent expansion of each ipl(t) is given by 

j oo 

= (t-m) k i7r)2 ^' 

1 ' 3=0 



where ipkj are the coefficients of the Laurent expansion of ip\{t). Next we 
substitute r = ^j, where we re-expand log A = 5^ >0 A;<5 2 ', this gives 



1 oo 



Inserting this expression into (I4.33I) we note that S k and S k cancel each other 
and we obtain an expansion in even powers of S 



n—1 n—1 



E sk tiw = E - k E ^'(E ^ 2l ) 2] ~ k ^r 



k=2 k=2 j=0 l>0 
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Since X(5) is even, analytic at 5 = and A(0) = 1 we use that X m expands 



DC 



fit 



A' 



to obtain 



=0 / 1=0 

n— 1 _ oo oo 



2/ 



n— 1 n— 1 oo oo 



fc=2 fc=2 j=0 «=0 

We will now change the summation indices in two steps to simplify this 
expression. First we let m = j + I which allow us to write the sums as 



n-l 



n-l 



E^) = EE* 2m E 

fc=2 m>0 j+l=m 



k=2 



Next we let p = 2j — which gives us the desired result 



n-l 



E^w = E* 2m E E ^ A 



p.i 



k=2 



m>0 k=2 2j-k=p 

l+j=m 

2m-2 

E* 2m E ^$ P . 

m>0 P=— n+1 



The derivation of the Laurent expansion of the second component of f)4.32p 
is analogous. The result is 



n-l 



[XI 



n-l 



/ 2m-2 \ 

m>0 p=— n+1 



fc=2 



2m-3 
\ m>0 p=— n+1 



(4.34) 



/ 



This finishes the derivation of the Laurent expansion of the formal separatrix 
around the singularity t = in. 



5 Close to identity maps 

In this section study an analytic family of close to identity maps of the form 

F 5 (x) = x + <5G 5 (x) (5.1) 



3 Do we need all these properties to draw this conclusion? 
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where 5 is a small parameter. Expanding a solution of the differential equa- 
tion x = Go( x ) into Taylor series in time, we can easily check that Fg is 
approximated by a map which shifts a point along a trajectory of the differ- 
ential equation by time 5: 

Fs = $ 5 Go +0(5 2 ). 

We say that Go generate the limit flow associated with the map Fg. 

If Fg is a family of area-preserving maps then the limit flow is divergence 
free and has a (possibly local) Hamiltonian function H , i.e., 

G = JV#o 

where J is the standard two dimensional symplectic matrix. 
5.1 Fixed points and their multipliers 

Theorem 5.1. If po is a hyperbolic saddle of the limit flow Gq then there 
is a unique smooth family pg of fixed points of Fg which converges to p as 
8 — > 0. Moreover, the fixed point depends analytically on 5 and is hyperbolic 
with multiplier \g. The multiplier depends analytically on 5 and 

\ 5 = l+5fi + O(5 2 ), 

where fio is an eigenvalue of -DGo(po) ■ 

Proof. The theorem easily follows from the implicit function theorem. In- 
deed, by the definition of a fixed point we have 

p 5 = F s (p s ), 

which is equivalent to 

G s (ps) = (5.2) 

due to the form of the map (15. II) . Since Fg is an analytic function, Gg is 
analytic. Since po is an equilibrium of the limit flow we have 

G (po) = 0. 

Since po is a saddle equilibrium of Go one of the eigenvalues of DGo(po) is 
positive and the other one is negative. Therefore 

detDGo(po) < 0, 
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and in particular det-DG (Po) 7^ 0. Then by the implicit function theorem 
there exists an analytic solution of (I5.2p for all 5 such that \S\ < 5 . 

Let us now consider the eigenvalues and eigenvectors of the fixed point. 

Let 

As :=DF s (ps) = id + 5DG s (p s ). 

The eigenvalues of -DGo(Po) are simple, hence eigenvalues and normalised 
eigenvectors of DGg(pg) are analytic in 5. The matrices DGg(pg) and Ag 
have the same eigenvectors. Moreover, if fj,(5) is an eigenvalue of DGg(pg), 
then Xg = 1 + 5p,(5) is an eigenvalue of Ag. □ 

5.2 Formal interpolation by a flow 

The results of this subsection are of independent interest. We show that the 
map Fg can be formally interpolated by an autonomous Hamiltonian flow. 

Proposition 5.2. Let Fg be an analytic family of area-preserving maps (15. ip 
defined on a simply connected domain D G M 2 (or D C C 2 ), then there is 
a formal Hamiltonian H$ with coefficients analytical on D, which formally 
interpolates F$. The Hamiltonian Hg is defined uniquely up to addition of a 
formal series in 5 only. 

Proof. For a function Fg let us denote by [Fg] n the n th -order coefficient of 
the Taylor series in S. In coordinates we write 

[Fg] n = (f n (x,y),g n (x,y)). 

We construct the formal series 

H 5 {x,y,5) = J2 6kH k(^,y), 

such that for every n 

[$* h = [F s )t for all 1 < i < n (5.3) 
where is the Hamiltonian flow generated by the Hamiltonian 

n-l 
k=0 

It is convenient to write the Hamiltonian flow generated by a Hamiltonian 
H'g using Lie series 

n=l 
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where is the Lie derivative generated by the Hamiltonian Hg. The Lie 
derivative is linear in H™. Then taking into account (15.41) we get 



E^j J2 L Hni ...L Hnm (x) 



m=l ni+ri2+. ..n m =t—m 
0<m<n— 1 



(5.5) 



We construct if„ inductively. Let n = 0. Then (15.31) takes the form 

o 

Taking into account H 5 . 5 1) we rewrite it in the form 

Lh (x) = G . 

Taking into account that in coordinates x = (x, y) and 

'dH dH 



dy ' dx 



we see that (15. 3p is satisfied for n = provided Hq is the Hamiltonian of the 
limit flow. 

We continue by induction. Suppose H , . . . , H n _i are chosen in such a 
way that the condition (15. 3p is satisfied for some n > 0. We show that H n 
can be chosen in such a way that (15.31) is satisfied for n replaced by n + 1. 

We note that (15.51) implies that 



for 1 < £ < n 



and 



n+l 



n+l 



Wx). 



Then the induction assumption implies that 



for 1 < £ < n . 



For f = n + 1 the equality is achieved provided H n satisfies the equation 



Wx) 



n+l 



n+l • 



(5.6) 



This equation has a solution in T> provided the right hand side is divergence 
free. The solution is defined uniquely up to addition of a constant. 
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In order to check that the right hand side has indeed zero divergence, we 
make a simple observation. Let F$ and F$ be two analytic families of close- 
to-identity area preserving maps such that [Fs]e = [F$]g for all < i < n. 
Then 

div[F 5 - F 5 } n+1 = 0. (5.7) 
Indeed, let us write the maps in the coordinates: 

F 5 (x, y) = U + S k a k (x, y),y + J^ s %( x , v) I 

V fe>i k>i J 

and a similar formula for F$. The preservation of the area is equivalent to 
the system 

d x a x + dyh = 0, 

n 

d x a n+1 + d y b n+1 = -^2{ai,b n+1 „i}, n > 1. 

i=i 

The system was obtained by expanding the Jacobian of F$ in power series of 
S and comparing the result with 1. Since = ag and bg = be for 1 < i < n, 
we immediately see that 

dx^n+i + dybn+i = d x a n+ \ + d y b n+ i (5-8) 

which is equivalent to (15.71) . 

Setting F$ = $frn and applying the last statement we see that the right 

5 

hand side of (15.61) has zero divergence and consequently H n is defined by the 
equation. The theorem follows by induction. □ 



5.3 Approximation of the local separatrix 

In this section we prove that the local separatrix of the map F$ is close to 
a local separatrix separatrix of any Hamiltonian flow which is 0(5 n+1 ) close 
to F$. The approximation is of order 0(5 n ). We note that the existence of 
a local separatrix near a hyperbolic fixed point of the map and a hyperbolic 
equilibrium of the flow follows from the Hartman-Grobman theorem^ 

The proof consists of two steps. First we transform the fixed point of the 
map to the origin and diagonalise its linear part by a linear symplectic change 
of variables. We also transform the Hamiltonian flow into a similar form. 

4 The proof does not really uses the fact that F is area-preserving and the flow is 
Hamiltonian. 
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Then we prove the theorem for the transformed map and the transformed 
flow. Then we come back to the original variables. 

Theorem 15.11 implies the map F§ has a fixed point ps which is analytic in 
5. Consider the following change of variables 

where B,5 is a symplectic matrix which diagonalises DF$(ps). This change of 
variables transforms F$ to 

F s = S- 1 oF 5 oS. 

The change of variables is chosen to ensure that the fixed point is moved to 
the origin 

F 5 (0) = 0, 

and the Jacobian of the map at the origin is diagonal: 

DF <5 (0) = A (5 :=diag(A ( 5,A,- 1 ). 

Note that the multipliers of the fixed point are left unchanged under any 
smooth change of variables. After the change of coordinates the unstable 
separatrix is parametrised by the function 

iP~ = So $~ 

which satisfies the equation 

$J{t + log \s) = F s $J{t)). (5.9) 

We note that the vector field generated by H$ has an equilibrium which is 
close to p<5 but does not necessary coincides with it. In general the linear part 
of the flow is not exactly diagonal after the substitution S. So we introduce 
another change of variables 

S:xH qj + Bax, 

where B<5 is a symplectic matrix that diagonalises JV 2 H$(qs). Using the 
implicit function theorem and smooth dependence of eigenvectors in a way 
similar to Section 15.11 we conclude that 

ps = q 5 + 0(cT), (5.10) 
B 5 = B s + 0(6 n ). 
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Then we define a new Hamiltonian function 

Let (fig denote the separatrix solution of the differential equation 

^#? = JVtf? . (5.1i; 

The original separatrix solution can be easily restored by the formula 



Let us write $s '■= ^ >1 °-i A ^„ to shorten the notation. By the definition of the 
Hamiltonian flow, 

£?(f + logA,) = (5.12) 

It is easy to see that 

where we have taken into account that //"^logAj = 1 + 0(5 n ) implies 
We have got 

$ s = F s + 0(8 n+1 ), (5.13) 
which will play a crucial role in the proof. 
Now we define 

£(f):=V>7(*) "¥£(*)■ (5-14) 
It follows from the definition of £(i) and equations f)5.12p and (15.91) that 

£(t + log X s ) = F 5 {^{t) + £(t)) - h(ft(t)). (5.15) 

We are proving that there exists a £ (i) such that 

£(t) = 0(5"e 2 *). (5.16) 

First we rewrite the equation for £ in the form 

£(t + log As) - A,£(*) = R(t,m) , (5-17) 

where 

R(t, £(t)) := F s (ffi(t) + £(*)) - - Ast(t). 
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Expanding Fs in Taylor series around (pg(t) we get 

R(t,£(t)) = R Q (t) + R 1 (t)at) + R 2 m)), 

where 



(5.18) 



Ro(t) ■ =F 6 (ffi(t))-Q 6 (m))> 

R 1 (t) : = DF s (^(t)) - A 5 , (5.19) 

mm ■■ =o(\m)\n (5.20) 

where the upper bound for R 2 (£(t)) is a standard estimate for the remainder 
of the Taylor expansion. Next we show that 

Ro(t) = 0(5 n e 2t ) (5.21) 

and 

R 1 (t) = 0(6e t ). (5.22) 

We begin with (15.221) . We note that since Fs is close to identity we have 
Fs = id + 5Gs and consequently DF$ = id + 5DGs and in particular As = 
DF s (0) = id + 5DG S (0). Therefore 

DF 5 (x) -A s = 8 (dG s (x) - DG s (0)) = 0(5|x|) . 

Substituting (pg(t) instead of x and using the upper bound tpg(t) = 0(e*) 
as Re(t) — > —00, we see the upper bound (15.221) follows from the definition 
(EH- 

The derivation of the upper bound for Ro(t) is a bit more sophisticated. 
First we note that 

F s (0) = <£ 5 (0) = and DF 5 (0) = D$ 5 (0) = A 5 . 

Consequently 

F 5 (x)-$,(x) = 0(|x| 2 ). 

Substituting (pg(t) instead of x and taking into account the definition of R 2 
we conclude that the components of are bounded analytic function on 

Re(t) < ro- Moreover ipg is 2ni periodic and analytic, therefore each compo- 
nent of the vector Ro(t) is 2ni is also periodic and analytic in t. Consequently 
the maximum of is achieved at the line Re(t) = r IH Using (15.131) we 
get for Re t < r 



\e- 2t R {t)\ < sup 

Re(t)=r 



F S (m)) - Mm)) 



o2t\ 



<c 2 5' 



n+l 



5 The 2ni periodicity allows a substitution z 
modulus principle to the result. 



e* and an application of a maximum 
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which is equivalent to the desired estimate (15.211) . 

Now we observe that if the following sum converges 

e(t) = ^A^ 1 /(t-HogA 5 ), 

k>l 

then it solves the difference equation 

£(t + logA 5 )-A 5 £(t) = /(i). 
Indeed, substituting the series into the left hand side of the equation we get 
£(t + logA a )-A 5 £ = ^A 5 fc -7(t-(^-i)logA 5 ) 

k>l 

-J2^ k sf(t-klog\ s ) = f(t). 

k>l 

Thus instead of solving (15.171) we solve the "integral equation" 

at) = J2 A l lR (t -klog\ s ,at -klogXs)), (5.23) 
k=l 

where R(t, is given by (I5.18p . We will show that we can solve (I5.23f) 
by applying contraction mapping arguments in a ball in the Banach space 
of functions which are analytic, 2rci periodic in a half plane Re(t) < r and 
have bounded norm 

||e||= sup |e- 2 ^(t)| 

Re(i)<r 

where | • | stands for the Euclidean norm in C 2 . Let 

B = { Ml^ll < K5 n } (5.24) 

where K > and r G M. will be chosen later in the proof. We note that if 
£ G B then for all Re(t) < r 

\m\ < U\\e 2Ke{t) 

< K5 n e 2Ke{t) . 

We will show that there exist K and tq such that the ball B is invariant 
under the operator in the right-hand side of fl5.23f) and that the operator is a 
contraction operator on the ball B. Let us denote the operator in the right- 
hand side of (I5.23P by T . We begin by showing that the ball B is invariant 
under J 7 , i.e. J 7 ^) G B for any £ G B. Let £ G B then 

oo 

HO = E A " _1 i^ + R ^ + MO)\ t - klog x 6 - (5-25) 

k=l 
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Using the upper bounds for Ro(t), Ri(t) and i? 2 (£(t)) from equations (15.211) . 
f)5.22p and (15.201) respectively we obtain 



-^olt-fciog x s 
-^l£lt-fclog x s 
-^(Olt-fclogA* 



< Cl 5 n+1 e 2K ^Xj 2 \ 

< c 2 6e^Xf k U\\, 
= c 3 5 2n K 2 e ARc ^Xj 4k , 



where ci, c 2 , C3 are constants. Since the norm of A<5 does note exceed Xs, we 
can use the above estimates to get an upper bound for the left hand side of 
(E25D: 



11^(011 < c lC4 6 n e 2ro + c 2 c A 8 n e 6r °KX^ + C3C4A7 V^A^e 



n„3r 



-2r2n-l j/2> 



(5.26) 



where C4 = sup 



\s\<s 



is finite since 



l-Aa 



0(5 x ). To show that there 



exists a K, r such that ||.F(£)|| — f° r & U 8 < 5o we start by choosing 
K = 4ciC4e 2r °. Then by choosing r = |log^^- we obtain that the second 
term in the right-hand side is bounded by 

c 2 c 4 5 n e 3ro A"A- 1 < ^5 n K. 



Finally letting 5 = 
of (I5T26D 

c 3 c 4 Xj z 5 
for all 6 < 5q. These choices yield 



2 Ic^a' wee ^ na ^ ^ ne ^ as ^ term m the right-hand side 



4 



||^(0|| 

hence ) G 23. 

Next we show that T is a contraction. Let £i,£ 2 ^ <^> an d consider 



Then we have 



= \R x (t) - &(*)) + ^(t, 6(t)) - i? 2 (t, £ 2 (0)l 

< (c 2 5e Re « + c 3 ||i? 2 || c2 AcTe 2Re «) |£ x (t) - £ 2 (t)| . 



11^(6) " ^(6)11 < £ SU P - ^(^))l,_ fclogA5 e 2Re(t) 

Re(t)<r 



k=l 
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Taking into account the previously obtained upper bound for Ri and R 2 we 

get 

~ ^(&)|| < c 5 (c 2 e ro + c 3 \\R 2 \\ C2 Ke^S^Xf) \\^ - &|| , 

where c 2 , C3 and C5 are constants. By choosing r$ sufficiently small the con- 
traction constant 

a = c 5 (c 2 e ro + c 3 ||i? 2 || c2 Ke^d^Xj 1 ) < ]- . 

It now follows by the contraction mapping principle that there exists a unique 
solution £ G B which satisfies equation (I5.16p . 
Now we recall that 

tP~ = So 

and 

ip s =Soip s 
Taking the difference we get the estimate: 

1>l(t)-V$(t) = So^(t)-So^(t) = So^(t)-So^(t) + 0(6 n ) 
= S fe(t) -m))+ 0(5") = 0(5 n ). 

We have proved that there is ro G R and a separatrix solution of (11 .5h such 
that 

Vr(t) = V?W+0(« B ) (5-27) 

for all t with Ret < r . 

5.4 Extension lemma 

In the previous section we have seen that the local unstable separatrix of the 
flow and the map approximates each other well for the half-plane Re(i) < r 
for some r G BL In this section we show that we can extend the domain in 
which the separatrices stay close to each other. 
Consider two close to identity maps 

F = id + 5G, (5.28) 
F = id + 5(5, 

The following lemma gives a closeness result for the difference of their first 
0(<5 -1 ) iterates. 
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Lemma 5.3 (Extension Lemma). Suppose two close to identity maps F$ and 
F$ are both defined on a 5-independent convex neighbourhood of a set V C C 2 
and on this neighbourhood F$ = F$ + 0(5 n+1 ). Then for any r > and any 
T > there are two constants 5q > and K such that for any two points 
Zq, Zq G D such that Zf. := F^(zq) G D for all k, < k < T8~ x , the inequality 
\z® — Zq I < rS n implies that 

\z h ~ Zk\ < K5 n for all < k < T5" 1 

where Zk = Fg(z ). 

To prove this lemma we will need the following Gronwall type inequality. 

Lemma 5.4. Let b > and (a n ) be an increasing sequence of positive num- 

n-l 

bers. If z n is a sequence such that z n ^ a n + b z k> then z n ^ a n (1 + b) n ~ . 

k=l 

Proof. By induction. The statement is true for n = 1. Assume that the 
statement is true for all n < p, i.e. z n ^ a n (1 + b) n ~ l . Then for n = p + 1 
we have 



Zp+i ^ a p+1 + b^Zk ^ a p+l + b^a k (1 + bf' 1 

k=l k=l 



_6T 
-6) 



= a p+1 (l + 6) p 

By induction the statement is true for all n G N. □ 



Proof of Lemma 15.31 Let us estimate the difference of the N th iterates of the 
two maps (15.281) . It is easy to see that 



N-l 

zn = z + 5 2j G(z k ) 

k=0 



N-l 



and z N = z + 5 G(5 fc ) . 



fc=0 



Let G — G = 5 n R, where i? is bounded due to the assumptions of the lemma. 
Taking the difference and applying the triangle inequality we get 

N-l 

\zn - z N \ = \z - z \ + 5 \G(z k ) - G(z k ) 

k=0 

N-l N-l 

< \z -z \+5Ys \G{h) ~ G{z k )\ + 5 n+1 \R{h) 

k=0 k=0 
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Since the domain of G is convex the map is Lipschitz 

\G{z k ) - G{z k )\ < L\z k - z k \ , 
where L = \\DG\\. Moreover R is bounded and we conclude 

jV-1 

\zn — z n I < 1 5 — zq\ + 5L \z k — z k \ + 5 n+1 N \\R\\ . 

k=0 

We recall that b := 1 5 — I < r $ n - The sequence ajy '■= S n r + 5 n+1 N\\R\\ 
is increasing and we apply the Gronwall type inequality of Lemma 15.41 to 
obtain 

|% - z N \ < (5 n r + 5 n+1 N \\R\\) (1 + SL)*' 1 . 
Then we use the textbook inequality 

(1 + 5L)^ < e 

as an upper bound for the last parenthesis in the right-hand side and get 

\zn - z N \ <5 n {r + 5N \\R\\) e"^" 1 ). 
Then for all < N < ? we have 

\zn -z n \ < 5 n (r + T|| J R||)e LT . 
Let K = (r + T\\R\\) exp(LT) to complete the proof. □ 

5.5 Application of the extension lemma 

We will now use Lemma 15.31 to extend the domain in which the local sepa- 
ratrix of the flow approximates the local separatrix of the map. We extend 
the approximation to the domain defined by 

T = {t E C : <p (t -s)GDVs>0, Re(t) < T} , (5.29) 

where <£> is the separatrix of the limit flow /xo0o = Go(cpo). Note that % is 
independent of 5. 

Corollary 5.5. Let ipj(t) and <p$(t) be defined as in Theorem \2.3\ and let 
T G M be fixed. Then ipj is defined on T and 

rl>i(t) = <ftS(t)+0(S») 

uniformly in %. 
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Remark 5.6. Uniformly in T means that 

\^(t)-^(t)\<K(D,T)6 n . 

Proof. By definition 

^(t + \ogX s ) = F s (4>J(t)) 

and 

^(t) = ^i JVfl ?(^(*))' 

which implies that 

(f + logA,) = *£&y(tf(f)) = $^ logW (^W)- 
Theorem 12.31 also states that 

F 5 = $ 1 H n + 0(5 n+1 ), 

o 

hence 

Let us denote this map by F§. Then we have 

^(t + log Xs) = F 5 (^(t)). 

It now follows that Fg and F$ are two close to identity maps such that 
F s = F s + 0(5 n+1 ). Define 

z k :=F^(^(t )) 

and 

4 ■= ft,)). 

Since Re(io) < ?~o it follows from Theorem 12.31 that 

<0(5 n ). 



To estimate 1^ (i) — <£$(t) we do 



Rc(f)-r 
log A 4 



iterations of the maps Fs, Fs 



hence the number of iterates is bounded by ^£°1 where L is a constant. 
Finally we note that by the definition of % it follows that z* k belongsto a 
5-neighbourhood of T> for all < k < L ( T ' ro ^ . The corollary now follows from 
Lemma 15.31 □ 
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